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ABSTRACT 

H I spatial power spectra (PS) were determined for a sample of 24 nearby dwarf irregular galax- 
ies selected from the LITTLE THINGS (Local Irregulars That Trace Luminosity Extremes - The 
Hi Nearby Galaxy Survey) sample. The two-dimensional (2D) power spectral indices asymptotically 
become a constant for each galaxy when a significant part of the line profile is integrated. For narrow 
channel maps, the PS become shallower as the channel width decreases, and this shallowing trend 
continues to our single channel maps. This implies that even the highest velocity resolution of 1.8 
km s _1 is not smaller than the thermal dispersion of the coolest, widespread Hi component. The 
one-dimensional PS of azimuthal profiles at different radii suggest that the shallower PS for narrower 
channel width is mainly contributed by the inner disks, which indicates that the inner disks have 
proportionally more cooler Hi than the outer disks. Galaxies with lower luminosity {Mb > —14.5 
mag) and star formation rate (SFR, log(SFR (Mq yr^ 1 )) < —2.1) tend to have steeper PS, which 
implies that the H I line-of-sight depths can be comparable with the radial length scales in low mass 
galaxies. A lack of a correlation between the inertial-range spectral indices and SFR surface density 
implies that either non-stellar power sources are playing a fundamental role in driving the interstellar 
medium (ISM) turbulent structure, or the nonlinear development of turbulent structures has little to 
do with the driving sources. 

Subject headings: galaxies: dwarf - galaxies: irregular - galaxies: ISM - ISM: structure - ISM: lines 
and bands - turbulence 



1. INTRODUCTION 

Pervasive interstellar medium (ISM) turbulence regu- 
lates star formation (SF) by creating local density en- 
hancements and countering gravitational collapse. Hier- 
archical SF (e.g. Efremov & Elmegreen 1998a; Gladwin 
et al. 1999; Zhang, Fall & Whitmore 2001) is a manifesta- 
tion of a fractal, turbulent ISM. Kinematic motions from 
turbulence prevent giant molecular clouds from collaps- 
ing to stars on the order of a free-fall timescale (Mac Low 
& Klessen 2004). Turbulent motion is also the dominant 
contributor of the total mid-plane pressure in the solar 
neighborhood (Boulares & Cox 1990; Jenkins & Tripp 
2001). In addition, the stellar initial mass function may 
be primarily shaped by turbulent fragmentation (Padoan 
& Nordlund 2002). Fast decay (~ a crossing time across 
the driving scale, Stone et al. 1998; Mac Low 1999) of 
the turbulent energy suggests a continuous driving mech- 
anism is in operation. It has been suggested that among 
the various power sources for turbulence, such as magne- 
torotational instabilities (MRI, e.g. Sellwood & Balbus 
1999), gravitational instabilities (e.g. Wada et al. 2002), 
thermal instabilities (e.g. Kritsuk & Norman 2002), and 
stellar energy (e.g. winds, supernovae), supernovae dom- 
inate the energy input to the ISM (e.g. Norman & Fer- 
rara 1996; Mac Low & Klessen 2004). Nevertheless, as 
pointed out by Elmegreen & Scalo (2004), gravitational 
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energy, which has an energy input rate an order of mag- 
nitude lower than that of supernovae, may have a higher 
efficiency for conversion into turbulence. Indeed, recent 
hydrodynamic galaxy simulations (Bournaud et al. 2010; 
Hopkins et al. 2012) found that ISM turbulence can be 
driven by gravitational processes on scales larger than or 
comparable to the Jeans length, with stellar feedback on 
smaller scales being essential also in maintaining a sta- 
ble ISM cloud structure. Observationally, it remains to 
be seen how ISM turbulence is related to SF. 

Fourier transform power spectra, which characterize 
the relative importance of structures at different scales, 
have been extensively used as a diagnostic for ISM struc- 
tures. A power-law behavior of the power spectra of H i- 
emission line intensities was found in the Milky Way 
(e.g. Crovisier & Dickey 1983; Green 1993; Dickey et 
al. 2001; Khalil et al. 2006), Small Magellanic Cloud 
(SMC, Stanimirovic et al. 2000), and Large Magellanic 
Cloud (LMC, Elmegreen et al. 2001). The power-law 
power spectrum is usually attributed to ISM turbulence. 
In unmagnetized, incompressible Kolmogorov turbulence 
(Kolmogorov 1941), kinetic energy is injected at large 
scale (the driving scale), and significant dissipation by 
viscosity occurs only at small scales (dissipation scale). 
The scale range in between the driving scale and dis- 
sipation scale, which has a power-law energy spectrum 
(E(k) ~ fc~ 5 / 3 ), is known as the inertial range. The ki- 
netic energy injected on large scales cascades to small 
scales without much loss in the inertial range. 

In incompressible, subsonic turbulence, density fluc- 
tuations passively follow the velocity field, which obeys 
Kolmogorov scaling (Lithwick & Goldreich 2001). In 
other words, density and velocity have similar Kol- 
mogorov power spectra P ~ fc~ n / 3 , where k is the 
three-dimensional (3D) wavevector. However, the tran- 
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sonic (e.g. warm neutral medium, WNM) or supersonic 
(e.g. cold neutral medium, CNM) nature of the ISM im- 
plies compressible turbulence. The ISM turbulence also 
involves self-gravity and magnetic fields. Numerical sim- 
ulations (Kim & Ryu 2005; Kowal, Lazarian, & Beres- 
nyak 2007; Gazol & Kim 2010) suggest that the spatial 
power spectrum is sensitive to the sonic Mach number 
and the Alfven Mach number, in the sense that higher 
sonic Mach number and Alfven Mach number (i.e. weaker 
magnetic forces) lead to steeper velocity power spectra 
and shallower density power spectra due to the formation 
of stronger shocks. Krumholz & McKee (2005) proposed 
that the sonic Mach number M is a factor in determining 
the star formation rate (SFR oc M -0 32 ). 

Hi is an important component of the ISM. Espe- 
cially in gas-rich dwarf irregular (dlrr) galaxies, neu- 
tral H I usually dominates the baryonic component (e.g. 
Zhang et al. 2012). Observationally, the intensity fluctu- 
ations of individual H I channel maps are caused by both 
3D real space projection and velocity mapping. Lazarian 
& Pogosyan (2000, hereafter LP00) found that the inten- 
sity fluctuations within thin velocity slices (less than the 
turbulent velocity dispersion at the studied scale) of the 
observed position-position- velocity (PPV) data cubes are 
generated or significantly affected by the turbulent veloc- 
ity field, whereas intensity fluctuations in thick velocity 
slices are dominantly caused by density fluctuations be- 
cause the velocity fluctuations are averaged out. Stan- 
imirovic & Lazarian (2001) applied this velocity channel 
analysis technique to the H i-emission line data of the 
SMC. They found that the power- law power spectra be- 
come steeper with increasing velocity slice width, from 
which the power spectral indices were derived for both 
the density and velocity fields. 

To investigate the relationship between the spatial 
power spectra and SF, we present a velocity channel anal- 
ysis for a subsample of LITTLE THINGS dlrr galaxies. 
The paper is structured as follows. In Section 2 we briefly 
describe the H i-emission line data used in this work. The 
power spectrum variations with channel width are pre- 
sented in Section 3. Section 4 gives a comparison between 
our image-domain-based and the visibility-based estima- 
tion of power spectra in the literature. In Section 5, we 
give the indication of the trend that power spectra be- 
come shallower as the channel width gets narrower. Sec- 
tion 6 explores various correlations between power spec- 
tral indices and SF-related quantities. Discussion about 
obtaining velocity spectral indices from our observations 
is given in Section 7. The summary follows in Section 8. 

2. SAMPLE AND DATA 

The 24 dlrr galaxies studied in this work (Table 1) were 
selected from the LITTLE THINGS sample (Hunter et 
al., in preparation) which were drawn from the large op- 
tical survey of dwarf galaxies by Hunter & Elmegreen 
(2006). The LITTLE THINGS is a large NRAO Very 
Large Array (VLAp*!) project, which was granted nearly 
376 hours of VLA time in the B-, C- and D-array config- 
urations to perform 21-cm Hi-emission line observations 
of a representative sample of 41 nearby (D < 10 Mpc) 

** The VLA is a facility of the National Radio Astronomy Obser- 
vatory (NRAO), itself a facility of the National Science Foundation 
operated under cooperative agreement by Associated Universities, 
Inc. 



dlrr galaxies. The sub-sample of 24 galaxies is chosen 
to be relatively face-on (inclination < 55°), and these 
galaxies cover a large range of galactic parameters, such 
as integrated luminosity (— 18 < Mb < —10 mag), cen- 
tral surface brightness (18.5 < (1q < 25.5 mag arcsec -2 ), 
SFR surface density (—4 < log(£sFR (M yr" 1 kpc" 2 )) 
< —1.3), and atomic gas richness (0.2 < M gas /M* < 26, 
see Zhang et al. 2012). Note that all but three (DDO 
50, NGC 3738, NGC 4214, Table 1) of our galaxies are 
fainter than the SMC (M B = -16.35 mag, Bekki & Sta- 
nimirovic 2009). 

To better handle extended emission, LITTLE THINGS 
adopted the multi-scale CLEAN (msclean) algorithm im- 
plemented in the Astronomical Image Processing System 
(AIPS). In particular, four different scale sizes, i.e. 0", 
15", 45" and 135", were chosen for mapping both small- 
and large-scale emission (Hunter et al., in preparation). 
For this work, we use the Hi data cubes in the image- 
domain created with the robust = 0.5 weighting scheme 
in the AIPS task IMAGR. The cubes are cleaned down to a 
flux level of 2 times the rms noise, determined in line- free 
channels. After the cleaning and mapping processes, we 
use the task BLANK to remove the noise. Briefly speaking, 
we first convolve the cube to a spatial resolution of 25", 
then in this low-resolution cube only regions containing 
emission (> 2 - 3 a) in at least three consecutive chan- 
nels are considered as areas of real emission and the rest 
is blanked (set to zero), and then we apply this Master 
Blanking cube to blank our full-resolution cube. 

The typical spatial resolution of the ROBUST maps is ~ 
7". Thirteen galaxies have a channel width of 1.3 km s _1 , 
with a real velocity resolution of 1.8 km s _1 . The other 
11 galaxies have a channel width of 2.6 km s _1 , with a 
velocity resolution of 2.6 km s _1 . Limited by the shortest 
baseline of the VLA's compact array configuration, our 
observations are blind to emission from structures with 
angular scales > ~ 15'. For our sample, except for IC 
10, all the galaxies have an Hi extent (the largest scale 
for power spectrum analysis) well below 15'. The VLA 
primary beam is 32' at 21 cm. The data cubes in the 
image-domain were corrected for primary beam attenua- 
tion, and then de-projected using geometrical parameters 
(galactic center, position angle and axial ratio) derived 
from ellipse fitting to the velocity-integrated Hi maps. 
An intrinsic axial ratio of 0.3 (Hodge & Hitchcock 1966) 
was adopted when converting the axial ratios to inclina- 
tion angles. We point out that adopting a slightly differ- 
ent intrinsic axial ratio introduces negligible uncertain- 
ties to our results in this work. 

3. POWER SPECTRUM VARIATIONS WITH CHANNEL 
WIDTH 

We determined the two-dimensional (2D) spatial power 
spectra of velocity slices with different width, by gradu- 
ally rebinning the channels with bin size increasing from 
the single channel width of 1.3 km s _1 (or 2.6 km s _1 ) 
to a sum over whole line profiles. Line-emission channels 
with at least 20 per cent of the peak flux of the global 
H I velocity profile were used in the analysis. The power 
spectra were obtained using the Fast Fourier Transform 
algorithm. Note that the power spectra presented in this 
paper are the average within each annulus (1 spatial 
frequency unit wide) in the 2D wavenumber k (k x ,k y ) 
space. The finite size of the synthesized beam results in 
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a pronounced decline at high spatial frequencies of the 
power spectrum (Figure [T]). Also, toward the lowest spa- 
tial frequencies of the power spectrum, some galaxies ex- 
hibit an obvious up-bending trend compared to higher 
spatial frequencies, which may be caused by some large- 
scale symmetric structures in the galaxies. Therefore, we 
fit each power spectrum with a power law (P oc k@) for 
linear scales from 1.5 times the beam size up to the point 
where the power spectrum starts bending upward. 

For a given slice width, the power spectral indices of in- 
dividual velocity slices with the same width are averaged 
to get the average spectral index < (3 > for that width. 
Most of our galaxies have non-zero inclinations, which 
means that an individual channel map shows only part 
of the galaxy, due to the galactic rotation. Thus galac- 
tic rotation may restrict the lowest spatial frequencies of 
the power spectra of narrow channel maps. If the fitted 
power spectral indices with and without the lowest one or 
two spatial frequencies are significantly different, we re- 
moved the lowest one or two spatial frequencies from the 
power spectrum fitting. As an example, Figure [2] shows 
some 2.6 km s^ 1 channel maps (black contours) of DDO 
133 overlaid on the integrated map (grey scale image). 
Our result suggests that, at a given channel width, differ- 
ent channel maps have about the same power spectrum, 
although they may show different parts of the galaxy. 

The variations of < j3 > with increasing slice width 
are shown in Figure [1] The error bars in Figure [T] are de- 
termined as the standard deviation of /3 divided by y/n, 
where n is the number of slices that were averaged to- 
gether. In Figure [TJ we also present the average power 
spectra of the 2.6 km s _1 channel maps (thin black solid 
lines) and the power spectra of the velocity-integrated 
maps (thick red solid lines). The vertical dashed lines in 
Figure [T] mark the linear scales of the synthesized beam. 
Table 2 lists the linear ranges used in the power-law fit- 
ting, the spectral indices of the velocity-integrated and 
the 1.3 km s _1 and 2.6 km s _1 channel maps, the ap- 
proximate slice width at which the power spectra start 
getting shallower for smaller slice widths, and the veloc- 
ity spectral indices, discussed below. 

As shown in Figure [U the power spectra are described 
very well by a single power law on large scales. Except 
for the two nearest galaxies (i.e. IC 1613, IC 10), the 
robust linear scales of the derived power spectra are > 
100 pc. DDO 187 and F564-V3 have the smallest range 
(less than a factor of 5) of linear scales used in the power 
spectrum fitting. As the slice width decreases, the power 
spectra of F564-V3 and DDO 101 get steeper first, then 
get shallower. Except for F564-V3 and DDO 101, the oth- 
ers exhibit gradually steeper power spectra as the veloc- 
ity slice gets thicker, and the spectral index approaches 
a constant as a significant part of the line profile is inte- 
grated together. 

The turbulent velocity field is correlated with spatial 
scales, in the sense that the velocity dispersion on larger 
scales is usually larger than that on smaller scales. In in- 
compressible turbulence, the turbulent density field pas- 
sively follows the velocity field. In a narrow channel map, 
small-scale structures from independent gas elements on 
the line of sight combine to give relatively more power on 
smaller scales. This results in a shallower power spectrum 
for a narrower channel map. However, if the spatially- 
uncorrelated thermal velocity is non-negligible in a given 



channel map, then small-scale structures with turbulent 
velocity dispersions smaller than the thermal speed will 
lose their distinction. If the thermal velocity is much 
larger than the turbulent velocity dispersion even on the 
largest scales, then the thermal velocity component will 
totally smooth out the spatial correlation present in the 
turbulent velocity field. Therefore, the shallowing trend 
for power spectra of narrower channel maps is a reflec- 
tion of the spatially-correlated turbulent velocity field, 
modulated by the thermal velocity. 

According to the theoretical calculations of LP00, the 
power-law indices of turbulent density and velocity power 
spectra can be obtained by changing the width of veloc- 
ity slices of the PPV data cube. If the effective velocity 
slice width, which is jointly determined by the instrumen- 
tal channel width and the thermal velocity (see below), 
is significantly smaller than the turbulent velocity dis- 
persion on the studied scales, then the velocity slice is 
regarded to be thin, and the power spectral indices re- 
main constant as the slice width further decreases. If the 
effective velocity slice width is significantly larger than 
the turbulent velocity dispersion on the studied scales, 
then the velocity slice is regarded to be thick, and the 
power spectral indices remain constant as the slice width 
further increases. The spectral indices of thick slices are 
determined solely by the turbulent density fluctuations. 
In the f3-Av plot of Figure [1] the asymptote on the side 
of thick velocity slices suggests that density fluctuations 
determine the intensity fluctuations there. The shallower 
power spectra of intensity fluctuations in narrower ve- 
locity slices implies more influence from velocity fluctu- 
ations. 

4. COMPARISON WITH A VISIBILITY-BASED 
ESTIMATION OF POWER SPECTRA 

Dutta et al. (2009b) have determined Hi power spectra 
of a sample of dwarf galaxies using a different method. 
They used data from the Giant Metrewave Radio Tele- 
scope (GMRT), and the range of baselines roughly corre- 
sponds to those of the VLA D and C array configurations. 
They determined the power spectra directly from the ob- 
served visibilities in the uv plane rather than producing 
maps first and then taking the Fourier Transform and 
determining the power spectra. This avoids an unneces- 
sary step (converting visibilities into maps) since the uv 
data are already in the Fourier domain. Here we prefer 
to identify real emission from the target by working in 
the image plane. Also, the correction for inclination and 
primary beam can only be properly done in the image 
plane. 

The Dutta et al. (2009b) study has two galaxies in 
common with our sample: DDO 210 and DDO 155. They 
found that DDO 210 has power spectral indices of —2.1 
± 0.6 and —2.3 ± 0.6 for velocity-integrated and sin- 
gle channel (1.7 km s _1 ) velocity slices, respectively, and 
DDO 155 has power spectral indices of —0.7 ± 0.3 and 
— 1.1 ± 0.4 for velocity-integrated and single channel ve- 
locity slices, respectively. These are compared to our val- 
ues of -3.66 ± 0.05 and -2.17 ± 0.09 for DDO 210, 
and -3.97 ± 0.15 and -2.59 ± 0.12 for DDO 155. Re- 
call that the single channel of 1.3 km s _1 is narrower 
than the real velocity resolution of ~ 1.8 km s _1 . The 
power spectra of DDO 155 are not power-law for linear 
scales above 1.5 times the beam size, and thus the quality 
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Figure 1. Power spectrum variations as a function of channel width. The average power spectra of 2.6 km s thick channel maps (thin 
black solid line), and power spectra of velocity-integrated maps (thick red solid line) are shown in the upper panel of each galaxy plot. 
The power-law fit to the power spectrum is overplotted as thin blue solid line. The power spectra have been arbitrarily shifted vertically. 
The vertical dashed lines mark the linear scales of the synthesized beam. Note that finite size of the synthesized beam causes the decline 
of power on linear scales smaller than the beam size. Variations of power-law spectral indices with channel width (in velocity units) are 
shown in the lower panel of each galaxy plot. 



of power-law fitting is very poor (Figure [lj . We notice 
that the power spectra of DDO f 55 determined by Dutta 
et al. (2009b) show similar behavior on large scales. So 
DDO 155 may not be a good case for comparison. For 
DDO 210, our single channel spectral index is consis- 
tent with that determined by Dutta et al. within the un- 
certainty of their measurement. Our velocity-integrated 
spectral index of DDO 210 is much steeper than that de- 
termined by Dutta et al. (2009b). Based on figure 8 of 
Dutta et al. (2009b), the multi-channel-integrated power 
spectrum of DDO 210 is a very good power law for lin- 
ear scales above ~ 0.18 kpc, below which the spectrum 
becomes very noisy and (thus) flat. Dutta et al. included 



the noisy, flat part of the power spectrum in their power- 
law fitting, which leads to a much less negative spectral 
index. 

The key difference between the image-domain-based 
and visibility-based power spectra lies in the different 
ways used to avoid noise bias. Noise present in the vis- 
ibilities or images shallows the derived power spectrum. 
To reduce the effect of noise in the uv plane, Dutta et al. 
determined the power spectrum by correlating the visi- 
bilities at slightly different baselines for which the noise 
is assumed to be uncorrelated. This method works on 
the assumption that the angular extent of the galaxy is 
much smaller than the primary beam of the telescope. 
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Fig. 1— Continued 



The compromise present in this method is that, on the 
one hand, sufficient visibility pairs are needed to get good 
statistics and on the other hand, the baseline differences 
should be as small as possible (< 1 /9 , where 9 is the an- 
gular extent of the target galaxy) to have strong enough 
correlations between visibilities of different baselines. 

In the image plane, we distinguish between real emis- 
sion and noise primarily in two steps. First, in the dirty 
map, only peaks 2 a above the noise level are regarded 
as signal and MSCLEAN-ed. Second, as described above, 
only regions with flux > 2 - 3 a in at least three con- 
secutive channels in the smoothed cubes are regarded as 
real emission, and the rest is blanked. We found that 
the power spectral indices of both DDO 210 and DDO 
155 would be about —1.5 regardless of the channel width 
if we do not apply blanking to the cubes. This implies 
that non-blanked cubes are predominantly affected by 
noise on small scales, msclean has been shown to work 



excellently in recovering extended structures (Cornwell 
2008; Rich et al. 2008). Compared to the classical CLEAN, 
msclean, with its scale-sensitive nature, can clean down 
to the noise level without leaving skewed noise on the 
residual map. Furthermore, the recovered fluxes from 
msclean (and blank) were found to be consistent with 
those from single-dish observations (Hunter et al., in 
preparation). Nevertheless, the noise present in the uv 
data can definitely affect the non-linear deconvolution 
involved in the image restoration process. 

Given these different approximations and assumptions 
in the two methods, it is not obvious to tell which method 
is better in determining power spectra. In this work, how- 
ever, we point out that, 1) by including the deep VLA 
B-array data, the outer part (longer baseline length) of 
the uv plane is better sampled in our data; 2) the two 
galaxies in common all have non-negligible inclinations. 
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5. EFFECTIVE VELOCITY SLICE WIDTH AND 
MULTI-PHASE NEUTRAL H I 

The distinction between the thin and thick regimes de- 
pends on a comparison between the squared turbulent 
velocity dispersion on the studied scale <J 2 urb , and the 
squared effective velocity slice width SV 2 (LPOO) 

°fur&, i < $V 2 , thick (l) 
°turb, i » 5V\ thin 

The thermal motion (vt), which is spatially incoherent, 
acts like a smoothing along the velocity dimension. Ther- 
mal broadening, together with the instrumental channel 
width, determine the effective velocity slice width. As- 
suming a uniform sensitivity across individual channels, 
the effective velocity slice width is given by 

SV ~2(Au 2 /6 + 24) 1 / 2 (2) 

where Av and vt are the individual channel width and 
typical thermal velocity, respectively (LPOO). It is the 
effective velocity slice width, rather than the channel 
width, that determines the effective slice thickness, and 
thus the variation of power spectral indices. 

The thermal velocity of the gas restricts the minimum 
effective velocity slice width that can be achieved. For an 
isothermal gas with temperature T, the effective veloc- 
ity slice width remains nearly constant when the chan- 
nel width_j3ecomes smaller than the thermal velocity vt 
(y/2kBT / rmr) , and thus the power spectral indices do 
not change. Empirically, the median channel width at 
the point where the power spectra start to get shal- 
lower for narrower channel width is ~ 15 km s _1 for 



our sample galaxies (Table 2). A two-phase (WNM and 
CNM) description of the neutral Hi was suggested by 
Field, Goldsmith, & Habing (1969). Wolfire et al. (2003) 
demonstrate that the interstellar H I at a temperature of 
~ 100 K and of ~ 8000 K can coexist in thermal equilib- 
rium. Therefore, if there is only thermally stable WNM, 
the power spectral indices would not change for channel 
widths narrower than ~ 6.5 km s" 1 (the thermal velocity 
of Hi gas with a temperature of 8000 K). For the thir- 
teen galaxies with a channel width of 1.3 km s _1 , the 
power spectra keep getting shallower for narrower chan- 
nel width, down to 1.3 km s . This suggests widespread 
H I with a temperature ~ 600 K (corresponds to a ther- 
mal velocity of 1.8 km s^ 1 ) or even lower, considering 
that the real velocity resolution is 1.8 km s _1 . In the 
Milky Way, Heiles & Troland (2003) found that more 
than 48% of the WNM, which accounts for ~ 60% of 
the total Hi, is in the thermally unstable regime with a 
broad temperature range from ~ 500 - 5000 K. A tem- 
perature of ~ 600 K lies in the broad boundary between 
the WNM and CNM. 

The typical velocity dispersion from Hi emission line 
observations of nearby star- forming galaxies is ~ 10 km 
s _1 (e.g. Cote, Carignan & Freeman 2000; Leroy et al. 
2008), which is the combination of turbulent and ther- 
mal velocity dispersions. Therefore, ~ 10 km s _1 can 
be regarded as an upper limit on the turbulent velocity 
on galactic scales. Also, since the maximum tempera- 
ture T max of neutral H I is ~ 10 4 K, the minimum possi- 
ble turbulent velocity dispersion Oturb should be ~ 4 km 
s _1 (\/10 2 — v\ ). For a given scale, the thin (thick) 
regime is reachecf if the squared effective slice width is 
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Figure 2. Every fifth 2.6 km s _1 channel map of DDO 133 used in our power spectrum analysis is shown as black contours. The grey 
scale image in each panel is the integrated intensity map. The synthesized beam is indicated at the left-bottom corner of the upper-left 
panel. The contour levels are 2, 10, and 20 times the rms noise. 



much smaller (larger) than the squared turbulent veloc- 
ity dispersion (Inequality 1). In the fi vs. Av plot, two 
asymptotes are expected for the two regimes of thin and 
thick effective velocity slices. As mentioned above, the 
Hi gas with a temperature of ~ 600 K has a thermal 
velocity vt ~ 1.8 km s _1 which is the highest spectral 
resolution for our data. With vt = 1-8 km s _1 and Av 
= 1.8 km s _1 , Equation 2 gives an effective velocity slice 
width of ~ 5 km s . Since we do not see the asymptotic 
behavior towards the single channel width in the /3 vs. 
Av plot (Figured]), the single channel, and thus the effec- 
tive velocity slice width of 5 km s" 1 does not reach the 
thin regime, at lease for the thermally unstable WNM 
components with a temperature < 600 K. Based on In- 
equality 1, the turbulent velocity dispersion aturb, i for 
the neutral Hi gas with a temperature < 600 K is not 
much higher (i.e. by less than a factor of 3) than ~ 5 km 
s , which is mildly supersonic. 

The line profiles of our galaxies cover a velocity range 
from ~ 20 to 130 km s _1 . The velocity-integrated slices 
lie in the thick slice regime, as implied by the asymptotes 
seen in Figure [TJ As mentioned above, our galaxies start 
reaching the thick regime at a channel width of ~ 15 km 
s _1 on average. With vt = 6.5 km s" 1 (%nm, stable = 
8000 K) and Av = 15 km s _1 , Equation 2 gives an effec- 
tive velocity slice width of 22 km s _1 . Therefore, to be 
in the thick regime defined in Inequality 1, the turbulent 
velocity dispersion of the thermally stable WNM should 
be much smaller (i.e. by at least a factor of 3) than 22 km 
s~ 1 . Therefore, the maximum Mach number for the ther- 



mally stable WNM components should be smaller than 

~ 1 {<Jturb, i/VTwNM.stable! °tW&, I < 22/3). 

The intensity fluctuations of the velocity-integrated 
maps are determined by density fluctuations including 
all possible phases of Hi. On the other hand, as ex- 
plained above, the shallowing trend for power spectra of 
channel maps narrower than 6.5 km s _1 is only caused 
by turbulent velocity fluctuations of cooler Hi, including 
the thermally unstable components and (possibly) some 
CNM. The large H I self- absorption survey by Gibson et 
al. (2005) found a smoothly distributed, albeit fluffy, cold 
phase of Hi in the Milky Way Galaxy, which is in line 
with our finding that Hi gas with a temperature of 600 
K or even lower is widespread in dlrr galaxies. 

Zhang et al. (2012) found that star-forming disks of 
most local dlrr galaxies have been shrinking at least dur- 
ing the past ~ Gyr. Cold atomic gas is the precursor to 
molecular cloud formation which shows an almost linear 
correlation with the SFR in nearby star-forming galaxies 
(e.g. Leroy et al. 2008). Possible radial variations of the 
f3-Av relation would reflect the relative spatial distribu- 
tions of H I gas with different temperatures. We took one- 
dimensional (ID) power spectra of azimuthal scans at 
different radii (e.g. Elmegreen et al. 2003) for our galax- 
ies. Figure [3] presents the results for three of our galaxies 
(DDO 70, DDO 47 and NGC 4214). The ID power spec- 
tra in Figure [3] are the average of adjacent azimuthal 
scans with a radial extent of two synthesized beams. We 
used the same geometric parameters derived above for 
the extraction of azimuthal profiles. The black symbols 
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Figure 3. Variations of ID spectral indices as a function of channel width for three galaxies. The red triangles denote the spectral indices 
of azimuthal scans at an inner radius, measured in units of V-band disk scalelength. The green upside down triangles denote the spectral 
indices of azimuthal scans at an outer radius. The 2D spectral indices presented in Figure [T] are shown as black circles. 



correspond to the 2D power spectra as shown in Figure 
[TJ The red and green symbols, respectively, denote the 
spectral indices at inner and outer radii, in units of the 
V-band disk scalelength i?^. In an isotropic field, the 
ID power spectra of azimuthal profiles are expected to 
be shallower than the 2D power spectra by ~ 1 because 
of the reduced dimension (see the power spectra of inner 
radii in Figure [3]). The azimuthal scans at the outer radii 
have systematically shallower power spectra than those 
at the inner radii, which may signify a gradual change 
from 3D to 2D geometry, possibly due to the longer az- 
imuthal scans at the outer radii. Obviously, the inner 
regions exhibit a much stronger shallowing trend toward 
narrower channel width. This suggests that the shallow- 
ing trend observed for the 2D power spectra, and thus the 
structure associated with the cooler Hi, is dominantly 
contributed by the inner, more actively star-forming re- 
gions. The more widespread cool H I in the inner disk may 
be caused by a higher mid-plane gas pressure (Wong & 
Blitz 2002; Blitz & Rosolowsky 2004) or a higher average 
gas volume density (Gao & Solomon 2004) there. Braun 
(1997) studied the neutral Hi properties of a sample of 
11 nearby spiral galaxies. He found positive radial gradi- 
ents of the H I kinematic temperature, and that there are 
more Hi components with narrow (< 6 km s^ 1 ) emission 
line profiles in the inner disks. 

6. POWER SPECTRAL INDEX VS. SF 

The ISM structure in dlrr galaxies can be significantly 
influenced by stellar feedback. Shells and holes of up to 
~ kpc scales are often seen in the Hi gas distribution 
of dwarf galaxies (e.g. Sargent et al. 1983; Puche et al. 
1992; Young & Lo 1997; Walter & Brinks 1999; Ott et 
al. 2001; Muller et al. 2003; Simpson, Hunter & Knezek 
2005; Cannon et al. 2011). It has been shown that stel- 
lar feedback can provide sufficient energy to produce the 
observed shells and holes (Warren et al. 2011). Stellar 
energy, such as winds, ionizing radiation, and super- 
nova explosions, dominates the energy injection into the 
ISM (Mac Low & Klessen 2004) on the corresponding 
scales, at least in star-forming galaxies. By studying the 
visibility-based power spectra of seven dwarf galaxies, 
Dutta et al. (2009b) found that the galaxies with higher 
SFR surface density tend to have steeper H I power spec- 
tra, implying more large-scale structures of H I for more 
intense SF. With a much larger sample, we can explore 
the possible relationship between the power spectra and 
SF. 

For most of the galaxies (Table 2), the minimum lin- 



ear scales covered by our power spectra are > 100 pc, 
which is probably comparable with the expected disk 
scale height of dlrr galaxies (e.g. van den Bergh 1988; 
Staveley-Smith et al. 1992; Carignan & Purton 1998; 
Elmegreen, Kim & Staveley-Smith 2001; Banerjee et al. 
2011). A break in the power- law power spectrum of the 
H I emission distribution has been observed for a few 
galaxies (LMC - Elmegreen et al. 2001; NGC 1058 - 
Dutta et al. 2009a) . The power spectral index for smaller 
scales (compared to the break scale) is steeper than for 
larger scales by ~ 1, which is found both in the LMC 
(Elmegreen et al. 2001) and simulations of gas-rich galax- 
ies (Bournaud et al. 2010). The flattening of the power 
spectrum on larger scales is understood as the transi- 
tion from isotropic 3D turbulence on smaller scales to 
anisotropic 2D turbulence on larger scales (Bournaud et 
al. 2010). The break is about the gas disk scale height 
(Elmegreen et al. 2001; Bournaud et al. 2010). 
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Figure 4. The synthesized beam size (pc) is plotted against the 
velocity-integrated power spectral indices determined by fitting a 
linear scale range from 1 to 1.5 times the beam size. 



As described in Section 3, we did power law fitting to 
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the power spectra for linear scales larger than 1.5 times 
the synthesized beam. This is because the finite size of 
the synthesized beam results in a pronounced decline at 
high spatial frequencies of the power spectra (Figure [T]). 
Nevertheless, we expect the power spectra at high spatial 
frequencies, although being noticeably influenced by the 
beam, still partially reflect the intrinsic power spectra of 
the galaxy at the resolution limit. Figure [4] presents the 
relation between the beam size (in units of pc) and the 
velocity-integrated power spectra fitted over linear scales 
from 1 to 1.5 times the beam size. Our galaxies fall into 
two groups in Figure IH separated around a beam size 
of ~ 200 - 300 pc, below which the galaxies have more 
negative indices. The two groups are also roughly divided 
by a spectral index of —3. This may reflect the effect of 
the disk thickness, because anisotropic 2D turbulence on 
large scales is expected to have much shallower power 
spectra than isotropic 3D turbulence on small scales. 

The SMC, which is like the most luminous galaxies in 
our sample in terms of stellar and gas masses, does not 
have a break in the power spectrum between linear scales 
of 30 pc and 4 kpc (Stanimirovic et al. 1999). This lack of 
a break indicates that the maximum transverse scale is 
comparable to or smaller than the H I line-of-sight depth. 
Perhaps the 3D behavior of the power spectrum of the 
SMC across a large range of linear scales is a result of 
the tidal interaction with the LMC during the past few 
Gyr. The numerical simulations by Bekki & Chiba (2007) 
found that almost 20% of the SMC's gas may have been 
accreted by the LMC during their recent interaction, and 
this may explain the origin of the LMC's intcrmcdiatc- 
age stellar populations with distinctively low metallicities 
(e.g. Geisler et al. 2003). So it is possible that the SMC's 
gas disk may have been stretched during the recent in- 
teraction, resulting in a thicker disk. For the two nearest 
galaxies in our sample, IC 10 and IC 1613, the power 
spectra also can be well fitted with a single power law 
between linear scales of ~ 50 pc and ~ 2 kpc. IC 1613 has 
a tidal index ~ 0.9 (Karachentsev et al. 2004), which in- 
dicates that this galaxy has been significantly influenced 
by the neighboring galaxies (e.g. M 33). Our Hi kinemat- 
ics analysis (Oh et al. in preparation) suggests that these 
two galaxies, and also the very low luminosity galaxies, 
such as DDO 210 and DDO 155, have ratios of the maxi- 
mum rotational velocity to velocity dispersion < 2, which 
indicates that the gas distribution may be a relatively 
thick disk or even triaxial ellipsoid. We note that, Roy- 
chowdhury et al. (2010, see also Sanchez- Janssen et al. 
2010) found a mean intrinsic axial ratio of ~ 0.6 for the 
Hi disks of dlrr galaxies with Mb > —14.5 mag, which 
suggests very thick gas disks in faint dlrr galaxies. 

Figures [5] and [6] plot the power spectral indices 
against SF-related quantities derived from the far- 
UV (FUV, SFR FU v): SF efficiency of the atomic 
gas (SFRpuv/ A/atom), SFR surface density normal- 
ized to th^ area within one U-band disk scalelength 
(SFR/ttR^ ), and SFR surface density within the ra- 
dius where the Hi column density falls off to 10 19 
cm- 2 (SFR/7ri?Hi 2 ), for 1.3 and 2.6 'km s^ 1 thick slices 
and for the velocity-integrated maps, respectively. For 
the formula used to derive SFRpuv from the FUV, the 
reader is referred to Hunter et al. (2010) and Zhang et al. 
(2012). The total atomic gas mass M at0 m (1.34x Mm to 
account for He) used here was collected from single-dish 



observations in the literature (see Hunter & Elmegreen 
2004 for the references). The open and filled symbols re- 
spectively denote the galaxies with spectral indices mea- 
sured at high spatial frequencies (1-1.5 times the beam) 
more negative and less negative than —3, i.e. the two 
groups mentioned above. The black circles and red trian- 
gles respectively denote the galaxies with Mb > —14.5 
mag and Mb < —14.5 mag. There is no strong correla- 
tion in these plots. However, Figure [5] suggests that all 
the galaxies with log(SFR) > —2.1 have narrow channel 
spectral indices around —2, whereas galaxies with lower 
SFR exhibit a much larger scatter toward more nega- 
tive spectral indices. The vertical dashed line in panel 
(a) of Figure [5] marks the division at log(SFR) = —2.1. 
Similarly, in Figure |6l the spectral indices of velocity- 
integrated maps exhibit a larger scatter (toward more 
negative indices) in galaxies with lower SFR. 

Given that the spectral indices show the bimodality 
with the SFR but not the SFR surface density, the bi- 
modality may just reflect the correlation between spec- 
tral indices and the mass (or size) of galaxies. In Figure 
El the spectral indices of 2.6 km s _1 channel maps and 
spectral indices of velocity-integrated maps are plotted 
against £?-band absolute magnitude Mb- There exists a 
similar bimodality as that present between power spec- 
tral indices and the SFR. The vertical dashed lines in 
Figure [7] mark the division of the bimodality at Mb = 
— 14.5 mag. Brighter galaxies tend to have flatter power 
spectra. This could be a result of disk thickness, in the 
sense that brighter galaxies tend to have a higher ra- 
tio of radial length to thickness, and so are more 2D on 
average. The three galaxies which have the most nega- 
tive spectral indices (for both narrow channel maps and 
velocity-integrated maps) are CVnldwA, DDO 70 and 
DDO 133. The negative tidal indices of these three galax- 
ies (Karachentsev et al. 2004) suggest that they do not 
have significant interactions with neighboring galaxies. 

This ratio of radial length to thickness also affects the 
regularity of galactic structure because the thickness is 
approximately the turbulent Jeans length, and therefore 
the typical size for stellar complexes. When the radial 
length is large compared to the thickness, there is room 
in the disk for a lot of star-forming clumps, giving the 
disk a somewhat uniform appearance. But when the ra- 
dial length is small compared to the thickness, there may 
be only a few star-forming complexes with large stochas- 
ticity, which can lead to an irregular appearance of the 
disk. Thus the correlation we find between the slope of 
the power spectrum and the galaxy absolute magnitude 
(or total SFR) may explain the findings of Lee et al. 
(2007), who explored the distribution of the Local Vol- 
ume galaxies in the Ms-Ha equivalent width plane. They 
found a bimodality around Mb ~ — 15 mag in the sense 
that galaxies with lower luminosity exhibit a larger (by a 
factor of 2) scatter in Ha equivalent width than galaxies 
with higher luminosity (— 19 < Mb < — 15 mag). 

The lack of a strong correlation between the inertial 
range spectral indices and SFR surface density (and also 
SF efficiency) may be unsurprising. In the classical Kol- 
mogorov turbulence, the energy transfer occurs locally 
(cascade) in the wavenumber space in the inertial range, 
and the driving only affects the energy input at the top 
of the cascade. Unlike this ideal Kolmogorov turbulence 
which has a single driving source on the large scale, the 
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Figure 5. Power spectral indices of 2.6 km s _1 channel maps are plotted against the SFR determined from the FUV (a), SF efficiency (b), 
SFR normalized by the area enclosed by the radius of one V-band disk scale length (c), and SFR, normalized by the area enclosed by the 
radius where the Hi column density falls off to 10 19 cm" 2 (d). The open and filled symbols respectively denote the galaxies with spectral 
indices measured at high spatial frequencies (linear scales between 1 and 1.5 times the beam size) more negative and less negative than —3. 
The black circles and red triangles respectively denote the galaxies with Mg > —14.5 mag and Mg < —14.5 mag. The horizontal dashed 
lines mark the expected ID (P ~ fc" 5 / 3 ) and 3D (P ~ fc" 11 / 3 ) power spectral indices of a Kolmogorov-type turbulence. The vertical 
dashed line in panel (a) marks the division (log(SFRpuv) = —2.1) of the bimodality discussed in the text. We notice that there is a slight 
trend for the galaxies with Mg < —14.5 mag in panel (d), in the sense that higher SFR surface density corresponds to shallower power 
spectra, but the scatter is too large to be sure. 



driving of turbulence in the ISM may span a wide range 
of scales. The stellar energy can drive turbulence from 
parsec-scale supernovae up to hundred-parsec-scale su- 
perbubbles (Mac Low & Klessen 2004). Similarly, grav- 
itational energy may feed the turbulence through sub- 
parsec gravitational collapse, up to galactic-scale gravi- 
tational instabilities (e.g. Wada et al. 2002; Elmegreen 
et al. 2003; Agertz et al. 2009) and galaxy interactions 
(e.g. Elmegreen & Scalo 2004), etc. Therefore, the sin- 
gle power-law behavior of Hi power spectra suggests 
that the turbulence may be driven (by whatever energy 



sources) over a wide range of physical scales, otherwise 
the power spectra should exhibit features, such as a flat- 
tening at low wavenumbers, at the primary driving scale 
(e.g. Nakamura & Li 2007; Padoan et al. 2009). Super- 
novae are the largest contributers of energy input to the 
ISM on scales comparable to or smaller than the disk 
thickness (Mac Low & Klessen 2004). However, it was 
suggested that supernova driven turbulence alone can- 
not explain the broad Hi emission lines, at least in the 
outer parts of disk galaxies (e.g. Dib, Bell & Burkert 
2006; Tamburro et al. 2009). Furthermore, simulations 
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Figure 6. Same as Figure[5]but for velocity-integrated maps. The slight correlation found in Figure 5d is also present here in panel (d). 



(e.g. Agertz et al. 2009; Bournaud et al. 2010) found 
that gravitational instabilities alone can reproduce the 
observed power-law power spectra, the stellar feedback 
does not significantly change the ISM statistical prop- 
erties established by gravitational instabilities. In addi- 
tion, the multi-phase nature of the ISM implies that ther- 
mal instabilities (e.g. Field 1965) may also be an impor- 
tant driving agent of turbulence (e.g. Hennebelle & Au- 
dit 2007; Gazol & Kim 2010). Thus turbulence may be 
driven by many sources and the properties of this turbu- 
lence on scales smaller than the sources may be indepen- 
dent of these sources, so a correlation between spectral 
index and SFR surface density might not be expected. 

The slope of the power spectrum may depend on the 
nature of the driving force (no matter whether it is stel- 
lar or non-stellar). Recent simulations done by Federrath, 



Klessen, & Schmidt (2009) suggest that, at a given Mach 
number, compressive forcing can lead to steeper density 
spectra than solenoidal forcing (rotational, incompress- 
ible). Kim & Ryu (2005), Kowal, Lazarian, & Bcrcsnyak 
(2007), and Gazol & Kim (2010) found a correlation 
between power spectral index and Mach number from 
magnetohydrodynamic (MHD) simulations, in the sense 
that higher sonic Mach number leads to shallower density 
power spectra. This correlation was interpreted as the re- 
sult of more small-scale density structures generated by 
stronger shocks in supersonic flows. In the present ob- 
servations, the lack of a correlation between the power 
spectral indices and the SFR surface density may im- 
ply that stellar feedback and other energy sources share 
similar characteristics of driving on average. Burkhart et 
al. (2010) also found little correlation between SF activ- 
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Figure 7. Power spectral indices of 2.6 km s _1 channel maps (left panel) and velocity-integrated maps (right panel) are plotted against 
B-band absolute magnitude Mg. The open and /li/eti symbols respectively denote the galaxies with spectral indices measured at high spatial 
frequencies (linear scales between 1 and 1.5 times the beam size) more negative and less negative than —3. The horizontal dashed lines 
mark the expected ID (P ~ k~ 5 / 3 ) and 3D (P ~ A;" 11 / 3 ) power spectral indices of a Kolmogorov-type turbulence. The vertical dashed 
lines mark the division of the bimodality (Mg = —14.5 mag) discussed in the text. 



ity and Mach number in a 2D map based on kurtosis of 
H I line profiles in the SMC. They found that regions with 
the highest sonic Mach number lie around the bar. Chap- 
pell & Scalo (2001) arrived at a similar conclusion using a 
multifractal spectrum analysis of low-mass star-forming 
cloud complexes: there is little correlation between the 
geometrical properties of the gas and the level of internal 
SF. 

As we discussed above, higher-luminosity galaxies {Mb 
< —14.5 mag) in our sample may have 2D turbulence 
on the studied scales and 3D turbulence on unresolved 
scales, whereas lower-luminosity galaxies may have only 
3D turbulence because the disks are thick relative to their 
radial length scales. Therefore, it may be desirable to ex- 
plore the relationship between SFR surface density and 
power spectral indices for higher-luminosity and lower- 
luminosity galaxies separately. Figures 5 and 6 suggest 
that our claim of a lack of a correlation between SFR sur- 
face density and spectral indices is valid for galaxies with 
both M B < -14.5 mag or M B > -14.5 mag. The lack of 
a correlation for the more luminous galaxies makes sense 
according to the above discussion if their power spectra 
are dominated by 2D turbulence and stellar feedback has 
little effect on scales larger than the disk thickness. The 
lack of a correlation for less luminous galaxies, which 
have relatively thicker disks than more luminous galax- 
ies, suggests that local SF does not strongly affect the 
spectral index of turbulence. 

7. THE TURBULENT VELOCITY FIELD 

Unlike the spectra that reflect density fluctuations, ve- 
locity spectra (specific kinetic energy spectra) are di- 
rectly related to the turbulent energy distribution across 
different scales. Intensity fluctuations within channel 
maps are contributed by both the density and veloc- 
ity fields. However, the relevant importance of density 



and velocity fields changes with the amplitude of the 
density fluctuations. For example, intensity fluctuations 
simply follow the velocity field if the density field is con- 
stant. Our thick velocity slices, especially the velocity- 
integrated ones, have intensity fluctuations dominated 
by variations in the density field. According to LP00, 
the spectral index of (2D) thick velocity slices /3thick (e.g. 
/^integrated listed in Table 2) is equal to the 3D density 
spectral index /?density, provided the maximum trans- 
verse scale is comparable to or smaller than the line-of- 
sight depth. In an isotropic 3D turbulent medium, if the 
density spectral index < —3, the velocity spectral index 
/^velocity is equal to — 2/3 t hin — 9, where /3thin is the spectral 
index of intensity fluctuations within thin velocity slices, 
/^velocity is equal to -3 - 2(f3 thin - /3 t hick) if the density 
spectral index > — 3. 

As discussed in Section 4, the velocity resolution of 2.6 
km s _1 does not yet reach the thin slice regime. Assum- 
ing the 2.6 km s -1 velocity slice is thin (i.e. /?thin = /?2.6) 
and the line-of-sight depth is comparable to the maxi- 
mum transverse scale we studied, we can obtain the ve- 
locity spectral indices /3 ve iocity which are listed in Table 
2. Sixteen galaxies have velocity power spectra steeper 
than that of Kolmogorov turbulence (for which /^velocity 
would be —11/3 for 3D turbulence). The real velocity 
spectral indices can be even steeper because we do not 
yet reach the "thin" regime. If the line-of-sight depth 
is smaller than the maximum transverse scale, the mea- 
sured /?thick and /3thin should be subtracted by <~ 1 to be 
used for the calculation of /^density and /3 ve iocity, and thus 
the resultant f3 dcnsity would be /3 t hick - 1, and /Velocity 
would be shallower than those listed in Table 2 by <~ 
2 when the /3density < —3, whereas /3 V eiocity determined 
for the galaxies with /3density > — 3 is not affected (see 
the relations in the preceding paragraph). We point out 
that, if the steepening factor of 2D turbulence relative to 



H I Power Spectra and the Turbulent ISM of Dwarf Irregular Galaxies 



13 



3D turbulence is not 1, then we would not know how to 
obtain ^density, and, in the case of /3 t hick < -3, /Velocity 
when the line-of-sight depth is smaller than the maxi- 
mum transverse scales. 

MHD isothermal simulations (e.g. Kritsuk et al. 2007; 
Federrath et al. 2009; Gazol & Kim 2010) suggest that 
fluids with higher Mach number (thus stronger compress- 
ibility) have steeper velocity power spectra (and shal- 
lower density power spectra). However, the isothermal 
assumption adopted in most simulations may not be valid 
for the multi-phase H I for which thermal instabilities can 
be an important driving mechanism of turbulence (Hcn- 
nebelle & Audit 2007; Gazol & Kim 2010). We emphasize 
that the multi-phase nature of the H I makes the veloc- 
ity spectral indices derived here questionable, because 
the density power spectra determined from thick veloc- 
ity slices may invoke density fields of all different phases 
of Hi, whereas the turbulent velocity fields as reflected 
by the shallower power spectra within narrower veloc- 
ity slices only include contributions from the thermally 
unstable WNM and (possibly) some CNM. 

8. SUMMARY 

We have studied the H I power spectral index variations 
with channel width for a sample of nearby dlrr galaxies. 
The majority of the 2D power spectra cover more than 
one decade of linear scales, from <~ hundred pc to several 
kpc. The main results are summarized as follows. 

1. The power spectral indices asymptotically become 
a constant for each galaxy when a significant part 
of the line profile is integrated, consistent with 
the theoretical calculations of LP00. This indi- 
cates that density fluctuations, including all possi- 
ble temperature components of Hi, determine the 
intensity fluctuations of our "thick" velocity slices. 

2. Starting at a channel width of ~ 15 km s _1 on av- 
erage for our sample, narrower channel maps have 
shallower power spectra. The shallowing trend, 
which is caused in part by turbulent velocity dis- 
persions of the thermally unstable WNM and possi- 
bly some CNM, continues down to the single chan- 
nel maps (1.3 km s _1 ). This continuation indi- 
cates that, first, even the highest velocity resolu- 
tion of 1.8 km s _1 is not smaller than the ther- 
mal dispersion of the coolest H I (< 600 K) which 
is widespread in our galaxies; if it were, then the 
spectral index would remain constant. Second, the 
turbulent velocity dispersion of the coolest Hi (< 
600 K) probed at our highest velocity resolution is 
not much larger than ~ 5 km s _1 , which means 
that the turbulence in this phase of Hi is mildly 
supersonic. 

3. Toward narrower channel width, the ID power 
spectra of azimuthal profiles at the inner radii have 
a stronger shallowing trend than those at the outer 
radii, which implies that the shallower power spec- 
tra for narrower channel maps are mainly con- 
tributed by the inner disks, and thus the inner, 
more actively star-forming regions have propor- 
tionally more cooler H I than the outer regions. 

4. The power spectra of IC 1613 and IC 10, which are 
the two nearest galaxies in our sample, can be well 



fitted with a single power law between linear scales 
of ~ 50 pc and ~ 2 kpc. This suggests that the 
H I line-of-sight depth may be comparable with the 
maximum transverse scales in these two galaxies. 

5. Our sample galaxies exhibit a bimodality in the 
spectral indices versus Mb (and also SFR) plane. 
The division of the bimodality is at Mb ~ —14.5 
mag and log(SFR (M Q yr" 1 )) 2.1. Galax- 
ies with higher luminosity and SFR tend to have 
shallower power spectra with a smaller scatter, 
whereas galaxies with lower luminosity and SFR 
exhibit a much larger scatter toward more nega- 
tive spectral indices. The bimodality may signify 
that higher-luminosity galaxies tend to have bigger 
gas disks compared to their thicknesses, whereas 
lower-luminosity galaxies may be better described 
as thick disks or even triaxial ellipsoids. 

6. The inertial range spectral indices of single chan- 
nel maps and velocity-integrated maps are not cor- 
related with the SFR surface density. This may 
imply that either stellar and non-stellar energy 
sources can excite turbulence with about the same 
power spectral index, or non-stellar energy sources 
are more important in driving ISM turbulence. 
The single power-law behavior of the power spec- 
tra indicates that the ISM turbulence may be 
driven, from whatever energy sources (stellar or 
non-stellar), over a wide range of physical scales, 
otherwise we should see features, such as a flatten- 
ing of power spectra at low wavenumbers, at the 
primary driving scale. 

The multi-phase (i.e. different temperature compo- 
nents) nature of galactic neutral H I means that the 
power spectra determined for different velocity slice 
widths trace different temperature components of Hi. 
Therefore, determining the turbulent velocity spectral 
indices may be difficult, and more theoretical work tak- 
ing into account the multi-phase nature of neutral H I is 
needed. 
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Table 1 

Galaxy Sample 



Galaxy 


Other Names 


D 


M B 


log(SFR) 


log(M atomic gas ) 


Rm 


R D 


Incl. 






(Mpc) 


(mag) 


(M yr- 1 ) 


(Mq) 


(kpc) 


(kpc) 


(degrees) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


DDO 210 . 


Aquarius Dwarf 


0.9 


-10.38 


-4.07 


6.52 


0.78 


0.17 


47 


M81dwA . 


KDG 052 


3.6 


-11.46 


-3.14 


7.26 


2.04 


0.26 





CVnldwA 


UGCA 292 


3.6 


-12.16 


-2.68 


7.81 


2.22 


0.57 


41 


DDO 155 . 


UGC 8091, GR 8 


2.2 


-12.25 


-2.71 


7.13 


1.35 


0.15 


48 


DDO 187 . 


UGC 9128 


2.2 


-12.38 


-3.14 


7.37 


1.23 


0.18 


32 


DDO 53 . . 


UGC 4459, VIIZw 238 


3.6 


-13.43 


-2.34 


8.39 


3.14 


0.72 


43 


DDO 75 . . 


UGCA 205, Sextans A 


1.3 


-13.72 


-2.18 


8.01 


3.09 


0.22 


34 


F564-V3 . 


LSBC D564-08 


8.7 


-13.67 


-3.05 


7.56 


3.37 


0.53 


53 


DDO 70 . . 


UGC 5373, Sextans B 


1.3 


-13.74 


-2.59 


7.72 


3.22 


0.48 


24 


NGC 4163 


UGC 7199 


2.9 


-13.95 


-1.73 


7.34 


0.60 


0.10 





IC 1613 . . 


UGC 668, DDO 8 


0.7 


-14.17 


-2.23 


7.66 


2.66 


0.58 


38 


Haro 29 . . 


UGCA 281, Mrk 209, IZw 36 


5.9 


-14.39 


-1.87 


8.01 


4.52 


0.29 


42 


DDO 46 . . 


UGC 3966 


6.1 


-14.36 


-2.05 


8.35 


4.76 


1.14 


29 


DDO 133 . 


UGC 7698 


3.5 


-14.36 


-2.13 


8.23 


3.82 


1.14 


43 


DDO 63 . . 


UGC 5139, Holmberg I 


3.9 


-14.58 


-2.09 


8.33 


4.22 


3.09 





DDO 87 . . 


UGC 5918, KDG 072, VIIZw 347 


7.7 


-14.52 


-2.15 


8.49 


8.51 


1.43 


28 


DDO 101 . 


UGC 6900 


6.4 


-14.40 


-2.68 


7.10 


2.36 


0.93 


49 


DDO 43 . . 


UGC 3860 


7.8 


-14.75 


-2.03 


8.40 


5.82 


0.61 


46 


DDO 52 . . 


UGC 4426 


10.3 


-15.05 


-2.03 


8.57 


6.99 


1.32 


53 


DDO 47 . . 


UGC 3974 


5.2 


-15.17 


-1.83 


8.73 


10.97 


1.36 


19 


IC 10 .... 


UGC 192 


0.7 


-15.75 


-1.73 a 


8.16 


4.00 


0.40 


41 


DDO 50 . . 


UGC 4305, Holmberg II 


3.4 


-16.39 


-1.18 


8.99 


8.67 


1.10 


47 


NGC 3738 


UGC 6565, Arp 234 


4.9 


-16.70 


-1.45 


8.32 


5.42 


0.78 


46 


NGC 4214 


UGC 7278 


3.0 


-17.26 


-1.03 


8.91 


8.62 


0.75 


26 



Note. — (1) Galaxy name. (2) The other commonly used names in the literature. (3) Distance from D. A. Hunter ct al. (2012, in preparation 
and references therein). (4) B-band absolute magnitude. (5) Logarithm of SFR derived from the FUV luminosity (Hunter ct al. 2010; Zhang et al. 
2012). (6) Logarithm of the atomic gas mass (1.34xMhi) is collected from single-dish observations in the literature (sec Hunter & Elmcgrccn 2004 
for the references). (7) The radius where the Hi column density falls off to 10 19 cm -2 . (8) Disk scale length measured on the V-band images. (9) 
The inclination angles derived by fitting the iso-intensity contours on the velocity- integrated H I maps. 
a No FUV observations, so the SFR is derived from Ha luminosity 
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Table 2 

Power Spectra 



Galaxy 


M i n - - n „ 
lv±111 scale 

(kpc) 


Max.. --I ~ 
(kpc) 


< Bi q > 
- Ml. 3 ^ 




< Bn « > 


a 


/-^integrated 


CT /3 


Ch cnan ge 

(km s ) 


Pvelocity 


( L ) 


( 2 ) 


(3) 




I 5 ) 


W 


{<) 


(°) 




(10) 




DDO 210 . 


0.11 


0.68 


-2.17 


0.09 


-2.47 


0.11 


-3.66 


0.05 


15 


-4.05 


M81DWA 


0.20 


1.35 


-1.33 


0.04 


-1.38 


0.04 


-1.98 


0.04 


27 


-4.19 


CVnldwA 


0.38 


2.37 


-3.09 


0.13 


-3.47 


0.14 


-4.33 


0.09 


23 


-2.07 


DDO 155 . 


0.27 


1.65 


-2.59 


0.12 


-2.86 


0.16 


-3.97 


0.15 


6 


-3.28 


DDO 187 . 


0.11 


0.43 


-1.43 


0.05 


-1.48 


0.06 


-1.98 


0.04 


9 


-4.00 


DDO 53 . . 


0.23 


2.23 






-1.71 


0.05 


-2.48 


0.04 


15 


-4.54 


DDO 75 . . 


0.09 


1.67 






-2.15 


0.04 


-2.80 


0.03 


23 


-4.31 


F564-V3 . 


1.31 


5.09 


-2.30 


0.08 


-2.46 


0.08 


-2.74 


0.07 


9 


-3.56 


DDO 70 . . 


0.14 


1.61 


-2.97 


0.06 


-3.29 


0.08 


-3.94 


0.07 


12 


-2.42 


NGC 4163 


0.16 


2.54 


-1.73 


0.04 


-2.06 


0.05 


-3.32 


0.04 


24 


-4.88 


IC 1613 . . 


0.05 


2.37 






-2.15 


0.07 


-2.85 


0.02 


26 


-4.38 


Haro 29 . . 


0.39 


2.07 


-1.80 


0.05 


-2.06 


0.06 


-2.87 


0.07 


17 


-4.62 


DDO 46 . . 


0.32 


2.29 


-1.59 


0.04 


-1.85 


0.04 


-2.75 


0.02 


14 


-4.80 


DDO 133 . 


0.43 


4.61 






-3.34 


0.07 


-3.94 


0.10 


13 


-2.31 


DDO 63 . . 


0.22 


2.05 






-1.44 


0.04 


-1.90 


0.03 


28 


-3.93 


DDO 87 . . 


0.48 


4.50 






-2.17 


0.05 


-2.90 


0.03 


15 


-4.46 


DDO 101 . 


0.59 


3.06 






-2.05 


0.08 


-2.00 


0.10 


10 


-2.88 


DDO 43 . . 


0.65 


5.87 






-2.13 


0.06 


-2.56 


0.06 


13 


-3.86 


DDO 52 . . 


0.85 


4.93 


-1.59 


0.05 


-1.93 


0.05 


-2.54 


0.05 


14 


-4.23 


DDO 47 . . 


0.42 


3.91 






-2.30 


0.05 


-3.14 


0.05 


21 


-4.39 


IC 10 .... 


0.04 


2.46 






-1.45 


0.04 


-1.59 


0.02 


13 


-3.28 


DDO 50 . . 


0.25 


2.72 






-1.63 


0.02 


-2.25 


0.02 


26 


-4.25 


NGC 3738 


0.32 


3.04 


-1.82 


0.04 


-2.13 


0.05 


-2.75 


0.04 


13 


-4.24 


NGC 4214 


0.18 


9.26 


-1.92 


0.03 


-2.10 


0.03 


-2.31 


0.05 


14 


-3.43 



Note. — (1) Galaxy name. (2) The minimum linear scale in the power-law fitting to the power spectra. Min sca i e is equal to 1.5 times the 
beam size. (3) The maximum linear scale in the power-law fitting to the power spectra. (4) The average 2D power spectral index of the 1.3 km 
s _1 velocity slices. (5) The uncertainty of the power spectral index of the 1.3 km s _1 velocity slices. (6) The average power spectral index of 
the 2.6 km s~ 1 velocity slices. (7) The uncertainty of the power spectral index of the 2.6 km s -1 velocity slices. (8) The power spectral index of 
the velocity-integrated map. (9) The uncertainty of the power spectral index of the velocity-integrated map. (10) The approximate channel width 
smaller than which the power spectra start to get shallower. (11) The velocity spectral index determined under the assumption that 1.3 km s _1 
velocity slices reach the thin regime and the line-of-sight depth is comparable to the maximum transverse scales. 



